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The interactions between the dielectric polarization and the fluctuations of the strain (stress) 
tensor in relaxor ferroelectrics are shown to give rise to the anisotropy of the anharmonic P^-term 
in the Landau-type free energy, however, the harmonic P^-term is still properly described by the 
rigid spherical random bond-random field model. These are the essential features of the compressible 
spherical dipolar glass model, which is used to calculate the singularities of the specific heat near 
field-induced critical points. The results agree with recent high-resolution calorimetric experiments 
in PMN [110]. 
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I. INTRODUCTION 

Relaxor ferroelectrics (relaxors) exhibit a variety of 
physical properties which are interesting for numerous 
practical applications, such as tunable capacitors, ul- 
trasonic transducers, actuators, and pyroelectric detec- 
tors Sometimes relaxors are regarded as a subgroup 
of incipient ferroelectrics in view of the fact that they 
do not possess a polarized long-range ordered phase in 
zero applied electric field. However, in contrast to nor- 
mal incipient ferroelectrics, relaxors undergo a freezing 
transition into a nonergodic glass-like phase below the 
so-called freezing temperature. If the relaxor is slowly 
cooled in a nonzero electric field it will pass through 
a sequence of quasi-stationary states. Thus, in order to 
stay close to thermal equilibrium, the experimental scale 
should increase steadily as the temperature is lowered. 
The corresponding E-T phase diagram is shown in Fig. 1 
for the case of PbMgi/gNba/sOs (PMN) in a field along 
the [110] direction 2]. The solid line in Fig. 1 separates 
the field-cooled dipolar glass phase from the field-induced 
long-range correlated ferroelectric phase. Similarly, the 
dotted line represents the boundary between the ergodic 
paraelectric phase and the frozen-in nonergodic dipolar 
glass phase. On approaching this line from the right, 
the longest dielectric relaxation time t{E, T) diverges 
according to the Vogel-Fulcher law [1, 3] , reflecting the 
random character of the relaxor state. 

Dielectric experiments in PMN [111] [1, Q have shown 
that there is no frequency dispersion of the dielectric sus- 
ceptibility in the region above the solid line, indicating 
that the relaxation times are finite and the system is er- 
godic. The transitions across the solid line, indicated 
by the arrows, are all first order and are characterized 
by a jump in the polarization P = P{E,T). As one 
moves towards higher temperatures, the size of the jump 
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FIG. 1: Phase diagram for a relaxor with 6 < 0. The solid 
line separates the dipolar glass (DG) phase from the field- 
induced ferroelectric (FE) phase. The arrows indicate the 
direction of the field-cooling (-heating) process. Dashed line: 
Supercritical regime. Dotted line: Freezing line separating the 
ergodic dipolar glass phase from the nonergodic one. Open 
circles: Data from Ref. 



becomes smaller and finally disappears at a liquid-vapor- 
type critical point Tcp, Ecp, where the transition is sec- 
ond order Beyond this point, the relaxor is in a su- 
percritical state characterized by a smooth evolution of 
P{E,T) and of the field-dependent dielectric susceptibil- 
ity x{E,T) = {dP/dE)T. The dashed line marks the 
positions of the maxima of xiE,T) (Widom line). 



II. POLARIZATION-STRESS COUPLING IN 
LANDAU FREE ENERGY 

It had been suggested earlier Q that when dealing with 
quasi-equilibrium states as in Fig. 1, the relaxor can be 
described in terms of a Landau theory based on the free 
energy density 
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For simplicity we are dealing with a scalar order param- 
eter P — P{E, T), corresponding to the polarization vec- 
tor along one of the symmetry directions in the crystal, 
i.e., [100], [110], or [111] in a system with average cubic 
symmetry. For an oblique direction of the electric field 
E the quartic term should be written as 



independent of the orientation and x"^"* < 0- Conse- 
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where the summation over all Cartesian indices is im- 
plied. Thus for a given symmetry direction, h will be a 
function of bijki^ which are components of a fourth rank 
tensor b. 

The first term _Fo in Eq. ([T|) contains the contribu- 
tion of all other degrees of freedom such as electrons, 
phonons, etc. The coefficient a = a{T) is related to the 
inverse quasistatic field-cooled susceptibility Xiy namely, 
a = (eoXi)^^. The susceptibility xi can be calculated 
from the static spherical random bond-random field (SR- 
BRF) model of relaxor ferroelectrics [13] and has the 
general form, 



Xi 
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which is well known from the theory of spin glasses and 
was found empirically to hold in the case of relaxors 
Here, Q is the Curie constant and Tq a measure of the 
average interaction between the elementary dipolar enti- 
ties in the system. In relaxors, these are known to be the 
polar nanoregions (PNRs), which are formed below the 
Burns temperature [ll]. Finally, q ~ q{T) is the dipolar 
glass order parameter, which is nonzero at all temper- 
atures due to the presence of quenched random electric 
fields [isj. In zero applied field, the order parameter q is 
determined by the real solution of the following algebraic 
equation 



ed by ti 
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The parameter J is defined in terms of the variance J"^ /N 
of the infinitely ranged random interactions of a spin- 
glass type, and A the variance of local random fields. 
For PMN, the estimated values are J/k ^ 217 K and 
= 0.001, whereas Tq in Eq. ^ is of the order 
kTo = Jo^ 0.9 J [13. 

The parameters 6, c, ... in Eq. ([T]) are related to the 
nonlinear susceptibilities xa, Xsv: which are defined as 
usual by the expansion P/eq = Ps+XiF+X3F^+X5P^ + 
■ ■ ■ . In relaxors, the spontaneous polarization vanishes, 
thus by definition Pg = 0, and the Landau coefficient b is 
given by 6 = — Xa/C^oXi)- It should be emphasized that 
in general both b and X3 j as well as higher order Landau 
coefficients in Eq. ^ , depend on the direction of the field 
E due to the anisotropy term ([2]). 

The SRBRF model was originally introduced for an 
ideal isotropic relaxor system [3] in a rigid environment. 
Thus the nonlinear susceptibility X3 derived from it is 



quently, br 



> 0, i.e., b. 



rigid 



is a positive scalar. Exper- 



iments on various relaxors have shown that X3 can either 
be positive or negative, depending on theparticular sys- 
tem studied and on the field orientation [y, [l^ . 

In order to derive a more general version of the model 
capable of reproducing the observed anisotropy of the 
coefficient 6, we introduce a coupling between the polar- 
ization P and the strain tensor Uij (or the internal stress 
tensor Xij). This suggests that we should consider the 
stress dependence of the Landau coefficients in Eq. ([1]). 
Focusing on the P^-term, we first introduce a generalized 
Landau coefficient Uki = eQ^{xi^)ki- Next, by expanding 
a{Xij, T)ki to linear order in Xij, we replace the P^-term 
in F by 
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The partial derivative is related to the electrostriction 
tensor Q^ju, namely, [ll] 
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By adding the elastic energy, the free energy F acquires 
an additional term, which can be written as 
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Here, C is the elastic constant tensor and (P^)ij = PiPj- 
Minimizing Fx with respect to Xij at constant tem- 
perature and field, we formally recover the free energy 
(H)), however, the quartic term is now replaced by the 
general expression 1^ with 

bijkl — brigid^ijSki + Bijkl, (8) 
where the fourth rank tensor B is given by 



B = 2Q C Q. 



(9) 



In relaxors, the electrostriction effect is usually large and 
the magnitude of the tensor components Bijki may ex- 
ceed the value of b^igid- Obviously, the sign of Bijki will 
in general depend on the balance between the individual 
components of Qijki and Cijki- Thus the resulting value 
of b for a symmetry direction can either be positive or 
negative. 

The generalized Landau free energy ([1]) with a(T) 
given by the SRBRF model and the P^-term having 
the form and with quartic coefficients (|S]) given by 
Eq. ([5]), will be referred to as the Compressible Spherical 
Dipolar Glass (CSDG) Model. 

We can evaluate the coefficients Bijki for the cases 
where values of Qijki and Cijki are explicitly known. 
In Table I, these are listed for the PMN crystal us- 
ing the Voigt notation, i.e., Pim = Pn, etc. The 
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TABLE I: Values of dj and Qij used to calculate Bij from Eq. Q and B'^' from Eq. (ITO)l . 
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"Reference [Igl; ''Reference [191]; "^References ''estimated |l7j] 



value of (544 can be estimated from the relation [17| 
Q44 ~ (Qii ~ '9i2)/2. From Eq. © we can then cal- 
culate the coefficients Bn, B12, and ^44 (see Table I). 
For a symmetry direction p, where p refers to [100], [110], 
or [111], the Landau coefficient h — fol^l in Ec^. ^ can be 
expressed in terms of hij as follows: 

611001 = 5u; (10a) 
6[iiol = i(6n+ 612 -f 2644); (10b) 
= i(6n+26i2-f4&44)- (10c) 

We can now write b^^ = b^igid + B^^\ where brigid > 0. 
Thus the rigid model always yields an isotropic, posi- 
tive contribution to &[p1, however, the new term is in 
general anisotropic. The corresponding values of S'^l for 
PMN are listed in Table I. While fil^ool > 0, we can see 
that B["o1 and Sli"! are both negative. Thus, 6[i"ol > 
and X3^°°' < 0. On the other hand, if \B^'p^\ > b„g^d for 
p = [110] and [111], the values of 6["0l and wih be 

negative, implying Xa^^^' > ^-nd Xa^^^' > 0, respectively. 
It is interesting to note that Tagantsev and Glazounov 
[H observed X3 > 0, but X3 < in PMN, in agree- 
ment with the above conclusions. Quasistatic measure- 
ments of polarization versus field along [111] and [100] 
directions also agree with these results ^] . As shown be- 
low, the sign of fol^' has important consequences for the 
existence of field-induced critical points for fields along 
the direction [p]. 

III. FIELD INDUCED CRITICAL POINTS 

The temperature and field dependence of the dielec- 
tric polarization during a field-cooled (or field-heated) 
quasi-stationary process is calculated by minimizing nu- 
merically the free energy ([1}. This procedure automally 
selects the correct solution of the minimization condition 
{dF/dP)E,T = 0. The parameter a = a{E,T) is cal- 
culated from the SRBRF model i8], while b and c are 
treated as free parameters. To ensure stability, we as- 
sume that c > and consider two cases, 6 > and 5 < 0. 
For 6 > 0, P{E,T) is found to increase monotonically 
with E and decrease with T, and no critical singularities 
of the susceptibility can be expected. For & < 0, how- 
ever, P{E) makes a discontinuos jump at some value of 
E at low temperatures. As T increases, the jump be- 
comes smaller and finally disappears at the critical point 




FIG. 2: Field dependence of P{E) for a relaxor with 6 < and 
several values of temperature T close to the critical tempera- 
ture Tcp, obtained by minimizing the free energy ([T]). Note 
that these calculations are only valid in the ergodic region 
above the freezing line shown in Fig. 1. 



EcPjTcp, where the slope of P{E) is infinite. This is 
illustrated in Fig. 2 for 6 = —0.2 and c — 0.08, corre- 
sponding to PMN [110]. A similar behavior of P{T) had 
been obtained earlier for PMN [111], where b = —1/3 
and c= \b\ m. 

The coordinates of the critical point are determined 
from the relations [1, 

a{Tcp)^—; EcP^—Pcp, (11) 
20c 25c 

where Pcp — ^—36/ (10c) is the polarization at the crit- 
ical point. The critical exponents at the field-induced 
critical point differ from the usual mean field exponents 
for ferroelectrics in zero field |8|]. 

The E-T phase diagram for PMN corresponding to 
cooling in a field along the [110] direction is plotted in 
Fig. 1 using the data points from Ref. [2j]. Similar phase 
diagrams were obtained earlier for PMN [100] and [111] 
[i|2l|. For i;||[100] no critical point was found, but the 
phase diagram for the [111] direction was shown to be 
analogous to the [110] case, in agreement with the above 
estimates for b^^^ . A similar conclusion had been reached 
earlier by Zhao et al. (2^ . 

The existence of field-induced phase transitions and 
critical points has recently been confirmed in PMN [110] 
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FIG. 3: (a) Calculated temperature dependence of the sucep- 
tibility x{E,T) for a set of field values E/Ecp, as indicated, 
(b) Same, but for the singular part of the specific heat from 
Eq. gni). 

by measuring the specific heat using high-resolution 
calorimetry !2]. The excess specific heat ACe{T), which 
is due to the contribution of the dipolar degrees of free- 
dom, namely, PNRs can be derived from the free energy 
Eq. ([1} by applying the thermodynamic relation for the 
entropy S — ~{dF/dT)E- We can write S — So -\- Sdip, 
where Sq — —{dFQ/dT)E and the dipolar part Sdip is 
defined as the contribution of all P-dependent terms in 
Eq. (HI, 

Sd^p = - QaiP' + ifoiP' + iciP*^ + • . • ^ , (12) 

where ai = da/dT, bi = dh/dT, etc., and the condition 
{dF/dP)E — has been applied jSj]. The dipolar excess 
specific heat capacity at constant field is given by IS.Ce — 
T{dSdip/dT)E, and at constant polarization similarly by 
ACp = TidSdip/ dT)p. These two quantities are related 
by the standard thermodynamic relation 

ACe = ACp+Tx{E,T)[{dE/dT)p]\ (13) 

where x{E,T) = {dP{E,T)/dE)T is the field-dependent 
susceptibility. The partial derivative in Eq. (|13p can be 
evaluated from the equation of state E = aP + hP^ -\- 
cP^ + ■■■, i.e., 

{dE/dT) p^aiP + hP'^ + ciP^ + ■■■ . (14) 

To calculate the expression on the right hand side we 
would, therefore, need to know the temperature depen- 
dence of the coefficients a, 6, c, etc. In practice, a{T) 



FIG. 4: (a) Experimental data showing the specific heat 
anomalies occurring in PMN [f fO] for four selected values of 
the electric field, obtained by high-resolution calorimetry 
(b) Calculated temperature dependence of the singular part of 
the excess specific heat C^"^ for a relaxor with 6 = —0.2 and 
c = 0.08 and the same field values as used in the experiment. 
The remaining parameters are taken from Ref. ^g]. 

is known from the SRBRF model through the relation 
— (soXi)^^, and ai{T) follows from the temperature 
derivative of Xi(r). On the other hand, b{T) and c{T) 
could, in principle, be estimated from the nonlinear sus- 
ceptibilities as given by the SRBRF model and the Q 
and C tensors. In the following we will simply assume 
that b and c are effectively constant in the temperature 
range of interest, and thus bi and ci will be neglected. 

The calculated temperature dependence of x{E, T) is 
shown in Fig. 3a for a set of values < -E < 2Ecp- Here 
we used the parameter values of J, Jp , and A determined 
earlier from the dielectric data [Tj], and the remaining 
parameters were chosen as 6 = —0.2 and c = 0.08. At 
E = 0, the zero-field cooled susceptibility xi is recov- 
ered. For < < Ecp, x{E,T) exhibits a jump at 
the first order transitions and diverges at the critical 
point. For E > Ecp, however, x{E,T) is character- 
ized by rounded maxima, in accordance with the smooth 
behavior of P{E, T) in the supercritical regime. 
The first term in Eq. (|13p is readily shown to be 

ACp = -T QaaP + h^P"" + \c2P'' + • • • ^ , (15) 

where 02 = dai/dT etc. Again, the corresponding deriva- 
tives of &, c, ... are not known, and we will neglect them. 
The quantity 02 can, however, be calculated from a(T). 



5 



It shows a sharp peak around the static "freezing" tem- 
perature Tf = (J^ + Ay/^/k, but is rather smah else- 
where. We may conclude that ACp will be nonsingular 
for all values of E, T, however, its precise behavior could 
only be determined if the values of b{T) and c(T) were 
known. 

The temperature dependence of the singular part of 
ACe{T) can be calculated from Eqs. ^ and (ITi|) . 
namely, 

ACr'-rx(£;,T)(aiF)2, (16) 

and is displayed in Fig. 3b for the same set of parame- 
ters as in Fig. 3a. The positions of the singularities of 
AC^™^(T) coincide with those of xiE,T), however, the 
direction of the jumps is reversed due to the last factor 
in Eq. (HH). 

The experimental data for ACe{T) in PMN [110] ob- 
tained by high-resolution calorimetry 0] are shown in 
Fig. 4a for four discrete values of the electric field E. For 
comparison, the theoretical prediction for the singular 
part of ACe{T), calculated from Eq. (IT51) at the same 
field values, is plotted in Fig. 4b. The experimental val- 
ues for the critical field and the critical temperature are 
^ven by Ecp = 8 kV/cm and Tcp ^ 240 K, respectively 

The predicted behavior of AC^'"^(T) qualitatively 
agrees with the experimental values of ACe{T), how- 
ever, the experimental anomalies appear to be broader 



than the calculated ones. There are several reasons for 
this broadening, for example, finite size effects, structural 
inhomogeneities, and slow relaxation. 

It should be stressed that at zero field {E = 0), no 
anomalies in ACe(T) were found in the entire tempera- 
ture range studied, in accordance with the CSDG model. 
This contrasts with the so-called random field scenario, 
according to which PMN [110] is assumed to undergo a 
ferroelectric phase transition at = 0. This point has 
been discussed in more detail in Ref. [2[. 



IV. CONCLUSIONS 

We have shown that the electrostrictive coupling be- 
tween the dielectric polarization P and the strain or 
stress tensor fluctuations in a relaxor ferroelectric gives 
rise to an anisotropy of the P^-term in the Landau free 
energy. For a given symmetry direction of the applied 
field the effective Landau coefficient b may become nega- 
tive, thus leading to the field-induced critical points. The 
compressible spherical dipolar glass model predicts singu- 
larities of the dipolar specific heat near critical points, in 
agreement with_ high-resolution calorimetry experiments 
in PMN [110] [2]. 
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